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ABSTRACT 

Strong toroidal magnetic fields generated in stellar collapse can generate magneto¬ 
centrifugal jets in analogy to those found in simulations of black hole accretion. 
Magneto-centrifugal jets may explain why all core collapse supernovae are found to 
be substantially asymmetric and predominantly bi-polar. We describe two phases: the 
initial LeBlanc-Wilson jet and a subsequent protopulsar or toroidal jet that propagates 
at about the core escape velocity. The prompt LeBlanc-Wilson jets will produce an 
excess of neutron-rich matter and hence cannot be the common origin of supernova 
explosions; similar, but less severe problems arise with the protopulsar jet that may be 
alleviated by partial evacuation along the axis by rotation. The jets will produce bow 
shocks that tend to expel matter, including iron and silicon, into equatorial tori. This 
may help to account for observations of the element distribution in Cas A. 

A magnetic “switch” mechanism may apply in rare instances (low density and large 
magnetic field), with subsequent increase in the speed and collimation of the toroidal 
jet. The conditions that turn the magnetic switch “on” would yield a jet that propa¬ 
gates rapidly and with small opening angle through the star, depositing relatively little 
momentum. The result could be enough infall to form a black hole. A third, highly 
relativistic jet from the rotating black hole could catch up to the protopulsar jet after 
it has emerged from the star. The interaction of these two jets plausibly could be the 
origin of the internal shocks thought to produce 7 -ray bursts and could explain the 
presence of iron lines in the afterglow. Recent estimates that typical 7 -ray burst energy 
is ~ 3 X 10®*^ erg imply either a very low efficiency for conversion of rotation into jets 
by the Blandford-Znajek mechanism, or a rather rapid turnoff of the jet process even 
though the black hole still rotates rapidly. 

Magnetars and “hypernovae” might arise in an intermediate parameter regime of 
energetic jets that yield larger magnetic fields and provide more energy than the routine 
case, but that are not so tightly collimated that they yield failed supernova. 

Subjeet headings: supernovae: general — pulsars: general — ISM: jets and outflows — 
gamma-rays; bursts 
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1. Introduction 

The problem of core-collapse has been with us for over 40 years (Hoyle &: Fowler 1960). Imme¬ 
diately after the discovery of pulsars, it was reasonable to explore the issue of whether or not the 
rotation and magnetic fields associated with pulsars could be a significant factor in the explosion 
mechanism (Ostriker & Gunn 1971; Bisnovatyi-Kogan 1971; Bisnovatyi-Kogan & Ruzmaikin 1976; 
Kundt 1976). With typical dipole fields of 10^^ G and rotation periods of several to several tens of 
milliseconds implying electrodynamic power of only ~ erg s“^, a strong robust explosion 

seemed unlikely. Several factors have led to a need to re-examine this conclusion. The principle 
one is the accumulating evidence that core collapse supernovae are distinctly asymmetric. Aside 
from its famous rings, HST observations of SN 1987A resolving the debris show that the ejecta 
are asymmetric with an axis that roughly aligns with the small axis of the rings (Pun et al. 2001, 
Wang et al. 2002). Chandra X-ray Observatory (CXO) observations of Gas A show that the jet 
and counter jet and associated structure are observable in the X-ray (Hughes et al. 2000; Hwang et 
al. 2000) as well as the optical (Fesen &: Gunderson 1996; and references therein). The most direct 
evidence bearing on this topic is from supernova spectropolarimetry which shows that substantial 
asymmetry is ubiquitous in core-collapse supernovae, and that a significant portion show strong 
evidence for a single, wavelength-independent axis of symmetry (Wang et al. 1996; Wang et al. 
2001). Many, even perhaps most, core-collapse supernovae are bi-polar (Wang et al. 2001, 2002). 

The strength of the asymmetry observed with polarimetry is higher (several %) in supernovae 
of Type Ib and Ic that represent exploding bare non-degenerate cores (Wang et al. 2001). The 
degree of asymmetry also rises as a function of time for Type H supernovae that have retained 
their hydrogen envelopes (from ^ 1% to ^ 1%) as the ejecta expand and one looks more deeply 
into the core material (Wang et al. 2001a; Leonard et al. 2000, 2001). Both of these trends 
suggest that it is the inner machine, the core collapse mechanism itself, that is responsible for the 
asymmetry. The observed polarization requires significant asymmetry; axis ratios exceeding 2 to 
1. To impose the observed strong asymmetry in the final homologous expansion it is plausible that 
an axial flow must be established and maintained for at least several dynamical time scales of the 
outer mantle/envelope. This is the operational definition of a jet. 

It also is significant that GXO observations have determined that jets are routinely associated 
with pulsars, not only in binary systems like SS 433, but also isolated objects like the Grab and 
Vela pulsars (Weisskopf et al. 2000; Helfand, Gotthelf, & Halpern 2001). One interpretation of 
these observations is that the present-day jets in these young objects are vestiges of a much more 
powerful MHD jet era that occurred in the first few seconds of the protopulsar phase when the 
compact objects were still inside their progenitor stars as they first attain nuclear densities. During 
this period the transient values of the magnetic field and rotation could have greatly exceeded those 
observed today, and indeed, as we suggest here, the reduction of the rotation rate and magnetic 
field to their present-day values could have been the process of energy release that powered the initial 
explosion. 
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Asymmetries associated with neutrino emission may produce dynamical asymmetries, but it 
is not clear that they can account for the polarization observations. Neutrino asymmetries may 
yield a short-lived, essentially impulsive effect (Shimizu, et al. 1994; Burrows & Hayes 1996; 
Fryer & Heger 2000; Lai et al. 2001). Expansion and transverse pressure gradients can wipe out 
transient asymmetries before homologous expansion is achieved (Chevalier & Soker 1989). “Finger” 
asymmetries might be preserved, but it is unclear that they can reproduce the common feature of a 
single symmetry axis that is substantially independent of space and time (Wang et al. 2001, 2002). 
Sufficient neutrino impulse might be delivered to the neutron star to yield a substantial runaway 
velocity (Burrows & Hayes 1996; Spruit & Phinney 1998), but it is difficult to see how this impulse 
can be communicated in a substantial and permanent way to the final ejecta trajectories. 

Highly-resolved, fully three dimensional, adaptive grid numerical calculations (Khokhlov 1998) 
have, however, established that non-relativistic axial jets of energy of order 10^^ ergs originating 
within the collapsed core can initiate a bi-polar asymmetric supernova explosion that is consistent 
with the spectropolarimetry (Khokhlov et al. 1999; Khokhlov & Hoflich 2001; Hoflich, Khokhlov & 
Wang 2001). Some imbalance in axial jets can also account for pulsar runaway velocities, specifically 
velocities that are parallel to the spin axis (Helfand, Gotthelf & Halpern 2001, and references 
therein). While a combination of neutrino-induced and jet-induced explosion may prove necessary 
for complete understanding of core-collapse explosions, jets as computed by Khokhlov et al. are 
sufficient. In this paper, we will ignore the possibility of neutrino-induced explosions as we focus 
on the possibilities of jet formation and jet-induced supernovae. For purposes of discussion, we will 
consider the sort of structure that forms a stalled shock after core bounce. 

Further work on jet-induced supernovae has been done by Khokhlov & Hoflich (2001) and 
Hoflich, Khokhlov & Wang (2001). These papers present well-resolved, three-dimensional jets that 
propagate from within the region of the original iron core until breakout from the star or stoppage 
in a red giant envelope in a single calculation that does not require artificial halting, rescaling, 
and restarting the jet. A bow shock forms at the head of the jet and spreads roughly cylindrically 
around each jet. The stellar matter is shocked by the bow shock and acts as a high-pressure confining 
medium by forming a cocoon around the jet. The jets become long bullets of high-density material 
moving through the background low-density material almost ballistically. Spreading is limited by 
a secondary shock that forms around each jet between the jet and the material already shocked 
by the bow shock. The laterally expanding bow shocks generated by the jets move towards the 
equator where they collide with each other. The result is that the material in the equatorial plane 
is compressed and accelerated more than material in other directions (excluding the jet material). 
The result is that heavy elements (e.g. O, Ca) are characteristically ejected in tori along the 
equator. Iron, silicon and other heavy elements in Cas A are distributed in this way (Hwang et al. 
2000), and there is some evidence for this distribution in SN 1987A (Wang et al. 2001b). 

The presence of a red giant envelope can stop the jet and dissipate the effects of the asymmetry 
in the outer envelope, but the asymmetry will still be manifest in the core. Radioactive matter 
ejected in the jets can alter the ionization structure and hence the shape of the photosphere of the 
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envelope even if the density structure is spherically symmetric (Hoflich, Khokhlov & and Wang 
2001). This will generate a finite polarization, even though the density distribution is spherical 
and the jets are stopped deep within the star and may account for the early polarization observed 
in Type II supernovae (Leonard et al. 2000; Wang et ah 2001). The jets are unstable to Kelvin- 
Helmholtz instabilities, but the growth time is long compared to the jet propagation time and no 
significant modulation of the jet is observed. 


Faster jets tend to make weaker explosions by propagating so rapidly and narrowly through the 
star that there is relatively little time to generate the transverse shocks that cause the explosion. 
The opening half angle of the jet is approximately (Wheeler et al. 2000), 


Venv 1 1 '^env.S ro '^env ,8 

~-~ 0.1 rad-^ ~ 5 - 

'IJbow Vbow,9 Vhow,9 


( 1 ) 


where Venv is the speed of sound of the envelope and Vbow is the speed of the bow shock, which is 
less than that of the inflow velocity of the jet material. As Vbow approaches the speed of light, the 
jet will be very narrow and affect a very small volume of the star as it propagates out. This feature 
may be relevant to the outcome of the explosion process and to the possibility of making 7 -ray 
bursts, as we will discuss in §4.3. Radiative transfer effects within optically thin portions of the jet 
may lead to significant instabilities in faster, lower density, jets (Hoflich, private communication). 
For related work in the context of black hole formation, see MacFadyen & Woosley 1999; Aloy, et 
al. 2000; MacFadyen, Woosley & Heger 2000; Zhang & Woosley 2001). 


The question then arises as to the mechanism of the production of the jets in routine core 
collapse events. A significant role for rotation and magnetic helds is an obvious candidate. Any 
mechanism that purports to account for routine pulsar formation may involve large transient mag¬ 
netic fields, but must ultimately be consistent with the distribution of deduced dipole strengths 
of “normal” pulsars of ~ lO^^-lO^^ G and with the distribution of initial rotation periods. There 
is growing evidence for a considerable spread in these quantities (e. g. Kaspi, et al. 2001 and 
references therein). In a minority of cases, the final dipole field might be consistent with a value of 
~ 10^^ G, yielding a “magnetar” (Duncan & Thompson 1992). Growing evidence for neutron stars 
of that field strength has been obtained by RXTE and other facilities (Kouveliotou et al. 1999; 
Ibrahim et al. 2001). The nature of the birth event of a magnetar is a separate problem that is 
significant in its own right. The ensemble of normal pulsar and magnetar births must also be con¬ 
sistent with observed nucleosynthetic abundances, a potentially crucial constraint on jet-induced 
supernova models that we discuss in §4. 


Possible physical mechanisms for inducing axial jets, asymmetric super novae, and related phe¬ 
nomena driven by magneto-rotational effects were considered by Wheeler et al. (2000). The means 
of amplifying magnetic fields by linear wrapping associated with differential rotation in the neu¬ 
tron star and possibly by a - H dynamos was discussed. Attention was focused on the effect of the 
resulting net dipole field on the torquing of the infalling plasma and newly formed neutron star 
and on the creation of strong Poynting flux (in the form, initially, of ultrarelativistic MHD waves) 
at the time-variable speed of light circle in analogy to pulsar radiation mechanisms, albeit buried 
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deeply in the core-collapse supernova ambience. Allusion was made to the existence and possible 
role of the primary toroidal field that is expected to form in routine collapse where protoneutron 
star rotation rates are not expected to trigger a - 17 dynamos and exponential field growth. 

In this paper, we will explore the capacity of the predominantly toroidal field generated near 
the core-mantle boundary to directly generate axial jets within core-collapse conditions by analogy 
with magneto-centrifugal models of jets in AGN (Koide, Shibata & Kudoh 1997; Meier et al. 1997; 
Romanova et al. 1998; Ustyugova et al. 1999; Meier, 1999; Koide, Meier, Shibata & Kudoh 2000; 
and references therein). The magneto-centrifugal models in the literature usually focus on black 
hole conditions where the field is anchored in the disk and at infinity. The physical picture discussed 
in this paper may prove especially robust since the field is anchored in substantial, corporeal objects: 
the outer, still collapsing layers of the progenitor star and the collapsed object — a neutron star. 
We find that the production of a strong toroidal field, substantially stronger than the 10^^ G dipole 
field of a pulsar, is nearly inevitable, and strong axial jets driven by that field equally so. The 
mechanisms described here also may prove to have a close astrophysical analog in the production of 
supersonic magneto-centrifugal jets in young stellar objects where also there is an outer envelope of 
infalling matter and an inner, differentially rotating object (Konigl & Pudritz 2000; Lery & Frank 
2000; Ouyed, Pudritz & Stone 1997). 

Our working hypothesis is thus that jets leading to asymmetric supernova explosions are as¬ 
sociated with routine core collapse and pulsar formation in roughly 90% of the cases. Perhaps 
10% of the core collapse events would be associated with magnetars. Magnetar formation might 
be associated with exceptionally strong explosions and asymmetries, but this is not necessarily the 
case. An even smaller fraction, perhaps 10“^ to 10“^ (Scalo & Wheeler 2002a) might involve black 
hole formation and the production of y-ray bursts. 

In §2 we outline the circumstances by which a newly formed neutron star could yield MHD 
jets. In §3, we discuss the basic physics of MHD jets and their formation. Section 4 explores 
the application of the physical principles of jet formation to supernova conditions, how magneto¬ 
centrifugal jets could be generated in core collapse, and some of their expected properties. We 
sketch a scenario in which extreme values of the rotation and magnetic field could lead to the 
failure of the initial supernova explosion with the subsequent collapse of the core to form a black 
hole and the attendant possibility to form a y-ray burst. In §5 we summarize our conclusions. 


2. The Production of Jets from Neutron Star Formation 

To understand routine, pulsar-forming, strongly asymmetric supernovae, there must be a 
deeper understanding of the generation of jets during the process of core collapse to form neu¬ 
tron stars. 

We identify three stages where magneto-rotational processes may come into play in the collapse 


process: 
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(1) The initial collapse or LeBlanc-Wilson phase, hereafter referred to as the “LeBlanc-Wilson jet” 
or “LW jet” (LeBlanc & Wilson 1970). In this phase, the wrapping of field lines of the original, 
pre-collapse core field within the nascent neutron star can give rise to a magneto-centrifugal jet 
or rising magnetic “bubble” that forms along the rotation axis at i? ~ 5 x 10 ^ cm and p 10 ® 
g cm“®. Parameters typical of modern pre-supernova calculations imply that the LeBlanc-Wilson 
stage is probably too weak to form a strong jet and explosion with observed supernova energies. 
Constraints from the production of r-process material also argue that this phase must not be robust 
(§4.1). 

(2) The protopulsar phase. There will be further wrapping of field lines by linear amplification 
of the contracting, cooling, de-leptonizing protoneutron star that is differentially rotating with 
respect to the infalling matter and the uncollapsed outer portions of the star. The toroidal fields 
can form within the newly formed, contracting, neutron star, at the shearing boundary between 
the neutron star and the infalling matter, and within the infalling matter. Between 10^^ and 10^^ 
ergs of rotational energy are available during the contracting deleptonization phase. A substantial 
portion of this energy will go into a toroidal field generated by differential rotation with a strength 
~ 10^^ - 10^® G. The jets are launched by toroidal magnetic “springs” that evolve into a steady- 
state, Blandford-Payne MHD wind (§3.1), the energy and momentum of which eject the stellar 
envelope. This is essentially the earliest phase in the birth of what will be the pulsar remnant, and 
is characterized by a period of rapid spin-down due to the production of the bipolar MHD outflow. 
This is a natural environment for a magneto-centrifugal jet and may be the routine mechanism 
for jet-induced, asymmetric supernova explosions. We will refer to this as the “protopulsar jet” or 
“slow toroidal jet” phase. 

(3) In relatively rare circumstances, depending on the distribution of rotation rates of the pre¬ 
supernova cores, the new neutron star may attain a rotation period of ~ 1 msec after deleptonization 
where the shear can compete with the convective motions and an a -D dynamo becomes possible. 
The resulting exponential growth of radial and toroidal fields generates fields of strength 10^® G or 
more. The rapid growth of the field coupled with the large rotational energy reservoir associated 
with millisecond rotation could lead to especially large explosions. A mix of a magneto-centrifugal 
jet and the emission of a strong Poynting flux at the light cylinder could lead to power of ~ 
1052 gj,g g-i extreme cases. It is possible that such cases have been witnessed as “hypernovae,” 
but strong asymmetries can affect estimates of energetics (Hoflich, Wheeler & Wang 1999). We will 
return to this topic in the conclusions. Another alternative is that very extreme conditions lead to 
especially fast, narrow jets that actually put less energy into the star (Khokhlov & Hoflich 2001). 
We will refer to these as “fast toroidal jets.” In this case, the MHD explosion mechanism may fail 
in the first instance and lead to the formation of a black hole. The formation of a rapidly-accreting 
black hole may, in turn, lead to the formation of an even faster, highly relativistic jet that interacts 
with the first, slower jet to produce a 7 -ray burst. We will explore this possibility in §4.3. 

Wheeler et al. (2000) discussed the possible physical processes, associated timescales, and 
energetics that could lead to the production of pulsars, jets, and asymmetric supernovae and to a 
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10^^ G magnetar in more extreme circumstances in the collapse of a massive stellar core. Wheeler 
et al. ( 2000 ) noted that when the neutron star contracts and speeds up two significant things may 
happen. One is that the rotational energy increases. The energy becomes significantly larger than 
that required to produce a supernova and, especially for recent estimates (Frail et al. 2001; Kumar 
& Panaitescu 2001), to drive a cosmic 7 -ray burst. In addition, the light cylinder may contract 
from a radius large compared to the Alfven radius to a radius comparable to that of the neutron 
star. This could disrupt the structure of an organized dipole field and promote the generation of 
intense MHD waves. The frequency of the MHD waves, ~ ~ 10^ Hz, would always be less 

than the plasma frequency, so these waves would be reflected internally and perhaps also channeled 
up the rotation axis. Collimation of the flow of energy may thus be expected. 

A jet emerging from deep in the collapsing iron core of a star will emerge from the helium 
core in 5 to 10 seconds. Wheeler et al. (2000; see also Nakamura 1999) argued that a jet launched 
after the deleptonization contraction of the neutron star could arrive at the surface of the helium 
star at about the same time. This argument is incorrect. The 5 to 10 seconds associated with 
the loss of neutrinos in the deleptonization phase represents the characteristic time scale for the 
luminosity, but the luminosity is enhanced at later phases by the contraction which results in a 
higher temperature at the neutrinosphere. The time scale for the contraction of the radius and 
hence the spin up of the neutron star is actually considerably shorter, of order one second (Wilson 
& Mayle 1993). This means that any stronger jet launched during the contraction will catch up to 
and dominate any magnetic “bubble” or jet launched promptly during the initial collapse before 
the first jet can emerge from the helium core. If this results in a supernova explosion, then any 
r-process material ejected, however slowly, from the neutron star in the first phase will be at least 
partially ejected in the explosion. Furthermore, if the collapse results in the formation of a black 
hole with the generation of a third, relativistic, jet, this jet will, in turn, catch up with any jet 
matter ejected earlier. The time when this happens will depend on the delay between the first 
iron-core collapse and the second protoneutron star collapse to form the black hole (§4.3). 


3. Magneto-centrifugal Jets 

3.1. The Blandford-Payne Mechanism: “Spring” and “Fling” 

The steady-state theory of magneto-rotationally-driven outflows was discussed by Blandford 
and Payne (1982) in the context of magnetized disks. Numerical treatments of such collimating 
winds have been performed by Ustyugova et al. (1995), Ouyed, Pudritz & Stone (1997), Krasnopol- 
sky, Li & Blandford (1999) and others. In the steady-state model, a magnetic field anchored at 
the base in a rotating object transfers energy and momentum to the plasma frozen onto the field 
lines. Differential rotation between the inner and outer regions of the outflowing plasma causes the 
field lines to be swept backward in an open helix; and the overall rotation of this structure pushes 
plasma upward parallel to the rotation axis and slowly compresses it toward that axis. Eventually 



this creates a collimated bi-polar outflow at large distances from the central object. This process is 
sometimes referred to as the “fling” model, since it employs centrifugal action to initially accelerate 
plasma along field lines. The forces that eventually complete the acceleration and collimation are 
magnetic pressure upward parallel to the axis and hoop stress. 

A similar, but time-dependent, process was investigated by Uchida & Shibata (1985) and more 
recently by Kudoh, Matsumoto & Shibata (1998), Uchida et ah (1999), and Koide et al. (2000). In 
this case an initially poloidal magnetic held {e.g., uniform or dipole) threads a differentially rotating 
medium, such as a thick disk or torus. As the evolution proceeds, the differential rotation coils 
the held in a barber-pole fashion, lifting and pinching material along the rotation axis. Depending 
on the strength of the held and the differential rotation, the early evolution of this system can be 
quite dynamic. The held will not have had time to expand to a broad open structure and approach 
the Blandford-Payne steady state. Instead, it may wind quickly to very large toroidal strengths 
and then expand vertically along the rotation axis, dynamically lifting and pinching plasma along 
the way. For this reason, this more dynamical process is sometimes referred to as the “spring” 
mechanism for jet outhow. Both spring and hing, however, use the same magnetic pressure and 
hoop stresses to perform the hnal acceleration and collimation. While both processes were hrst 
discussed in the context of accretion disks, they will have close analogies in collapsing supernova 
cores, especially if the core rotation results in signihcant core battening and an accretion-torus-like 
structure surrounding the core. 

For non-relativistic outhow, the power output of an MHD how is given by (Blandford & Payne 
1982) 

Lmhd = / 2, ( 2 ) 

where B, D, and R are, respectively, the magnetic held at the base of the jet, the jet angular 
velocity, and the size scale of the region where outhowing matter is injected along the held lines. 


3.2. The Magnetic Switch 


If, in the above cases, the magnetic held is sufficiently strong or the rotation sufficiently fast, 
then the power output can exceed a critical dynamical value, given by the ratio of the escape energy 
to the free-fall time. For the case of a constant density corona of radius R, this can be expressed 
as: 


L 


crit 


Eesc 

Vf 


mej 


( 3 ) 


where the escape energy is taken to be GMmej/R- Simulations of magnetized accretion disk 
coronae with a height comparable to the radial extent {niej ~ AnR^p; see Meier et al. 1997; Meier 
1999) have shown that under these conditions much of the MHD power may be converted quickly 
into kinetic energy, creating a fast, highly-collimated jet close to the injection region, and that this 
condition can persist well into the steady state. In the current situation, the mass to be ejected 
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that determine the escape energy is that in a cone of fractional solid angle /ah = Ahl/dyr, where 
An is the solid angle of the jets emerging from the poles of the neutron star. Assuming that the 
mass in this cone is concentrated near the lower portions where the density is p at radius ii, we 
take mej ~ 4/37rii^/An and hence; 


L 


crit 
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Vf 
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( 4 ) 


For Lmhd < Lcrit^ the outflow will be a classical Blandford-Payne-type outflow, collimating slowly 
with the jet thrust spread over a broad region. (See, e.g., the steady-state simulations of Krasnopol- 
sky et al. 1999). This criterion for a slow jet can be expressed in terms of a dimensionless parameter, 
n, as 


u = 


< 1 . 


3n \ 

For V > 1, the jet may have a much higher velocity (with Vjet ^ va > Vesc) 
thrust concentrated in a very narrow axial spindle (Meier et al. 1997). 


/ Lmhd \ 

V Rcrit J 


_ VA r-1/2 
“ „ JAn 

^esc 


( 5 ) 

early on, with the jet 


4. Magneto-Centrifugal Jets in Supernovae 

We will now explore the application of the “spring” and “fling” and fast and slow jets to the 
supernova collapse problem and consider in more detail the three stages of the collapse that may 
lead to MHD jets as outlined in §2. 


4.1. The LeBlanc-Wilson Jet 

The first type of MHD outflow we consider is that discovered by LeBlanc & Wilson (1970; 
hereafter LW), also treated later by Symbalisty (1984). LW studied the magneto-rotational collapse 
of a 7 Mq iron core by numerically solving the two-dimensional MHD equations coupled to the 
equation for neutrino transport. Their simulations showed that a “circulation vortex” formed along 
the rotation axis, at about the time of core bounce, amplified the magnetic held there to of order 
~ 10^^ Gauss, and created two oppositely-directed, magnetically-driven, high-density, supersonic 
jets of material emanating from inside the protoneutron star. In the calculation of LW, the jet 
formed at about 50 km, in that of Symbalisty, at about 30 km. LW estimated that their jet carried 
away ~ 10^^ g with ~ 1 — 2 x 10^^ erg in ~ 1 s. Khokhlov et al. (1999) assumed characteristics of 
the LeBlanc-Wilson jet for their proof of principle calculations. 

The calculation of LW was extended by Meier et al. (1976) to treat stellar cores with mass, 
central density, and rotational characteristics more commensurate with the results of modern ad¬ 
vanced stellar evolutionary calculations. One of their main constraints on the LW mechanism was 
that it not over-produce r-process material in the Galaxy, a limit of roughly 2 x 10“^ Mq of highly- 
neutronized ejecta per supernova. Meier et al. pointed out that this constraint was not compatible 
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with the notion that the LW mechanism provides the explosion energy for most supernovae: if it 
were to do so, then the r-process elements would be over-produced by a factor of ~ 1000. Therefore, 
either the LW jet provides the explosion energy (and most of the Galaxy’s r-process material) in 
only a very few rare supernovae, or LW jets or bubbles are quite common in most core-collapse 
supernovae, but eject only ~ 10“^ M 0 per event, and are not important in the supernova’s overall 
energetics because the progenitor did not have especially large initial rotation and magnetic field. 

In the latter, more plausible, case, Meier et al. (1976) showed that the LW events would proceed 
as “buoyancy instabilities,” with the magnetized material rising subsonically up the rotation axis of 
the protoneutron star and bursting into the mantle above. In the more modern language of MHD 
winds and jets, these LW events will result in Blandford-Payne-type MHD wind flows emanating 
from the upper interior of the protoneutron star. These MHD “fountains” will be initially sub-sonic, 
relative to the core sound speed, and rather confined by the high pressure, but when they burst into 
the mantle, they may be supersonic with respect to the sound speed there. In a completely-bare 
core collapse, with little or no mantle to disrupt these weak outflows, these events eventually could 
accelerate to supersonic speeds and form jets of r-process material, but in most supernovae the LW 
ejecta are likely to quickly mix with the polar regions of the infalling mantle material. 

Taking the collapsing protoneutron star to be about 1.5 Mq in mass, with initial and final 
radii of ~ 10® cm and ~ 5 x 10® cm, respectively, we estimate from Meier et al. (1976) that the 
magnetic field still will be rather high (few x 10^® Gauss) in the circulation-vortex region near the 
axis, but that the differential rotation rate there will be rather low, < 100 rad s“^. From the 
r-process constraints and equation (5) we find that the magnetic switch parameter v will be of 
order 0.03 or less (for /ao ~ 1), nearly independent of the exact values of the collapse parameter 
and the efficiency of conversion of rotational energy to magnetic energy. This result that the MHD 
outflow is well below the magnetic switch point is consistent with the LW fountain being an initially 
subsonic, Blandford-Payne-like, polar outflow into the mantle. 


4.2. The Slow Toroidal Jet 

4-2.1. Generation of the Toroidal Field 

Rotation alone can inhibit the contraction of the newly-formed neutron star (Symbalisty 1984; 
Monchmeyer, et al. 1991; Fryer & Heger 2000). However the addition of magnetic fields provides 
a source of dissipation that will allow further contraction, thus enhancing the rotational energy, 
the differential rotation, and ultimately the field energy. The presence of even a small field, there¬ 
fore, can significantly increase the possibility of creating strong axial jets. As discussed in the 
introduction, Wheeler et al. (2000) considered the linear wrapping phase of field amplification due 
to differential rotation within the neutron star (Meier, et al. 1976; Kluzniak & Ruderman 1998; 
Ruderman, Tao & Kluzniak 2000), but they did not emphasize the resulting toroidal field, only the 
dipole component of that field (taken to be ~ 1% of the toroidal field) and its possible effects in 
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concert with the Alfven surface and speed of light circle. Here we argue that the principle focus 
should, rather, be on the strong toroidal field itself. By analogy to the work on the creation of jets 
by MHD processes in accretion disks (§3) the toroidal field can act as a “spring” to launch jets. 
We expect this sort of jet to overtake and supplant the original LeBlanc-Wilson jet or bubble, as 
illustrated schematically in Figure 1. 


The original rotational energy of the protoneutron star (PNS) is approximately, 

rp 1 r o2 n T nSO f ^NS \ f PpNS \ ^ f PpNS \ ^ 

and the rotational energy after the deleptonization and contraction is about, 
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( 7 ) 


It is not clear how much of the latter energy can be dissipated in the magnetic field, but there is 
substantial energy for a supernova for any final rotation period of less than about 7 msec. Strictly 
speaking only the energy associated with the differential rotation can be tapped, but the neutron 
star is always in strong differential rotation compared to the progenitor star. There may be constant 
torques on the neutron star until the explosion succeeds or fails, so this issue is complex and involves 
the magnetic fields we will now explore. 


If the pre-collapse progenitor core has a field strength comparable to that of a magnetized white 
dwarf, ~ 10® G, then a field of r\j 10^^ G could arise during the collapse simply from flux-freezing. 
This field can be amplified further by differential rotation in the neutron star (Meier et al. 1976; 
Kluzniak & Ruderman 1998) to produce a strong toroidal field. For a collapsing white dwarf. The 
protoneutron star is likely to be differentially rotating with angular velocity H increasing outward 
since the protoneutron star will be “stiffer” than the highly degenerate white dwarf progenitor 
(Ruderman, Tao & Kluzniak 2000). For a neutron star formed by homologous iron core collapse, 
the angular velocity of the homologous core will increase outward within the homologous core at 
core bounce, but after bounce the angular velocity is nearly constant in the inner core and decrease 
with radius beyond the boundary of the initial homologous core (Akiyama, Wheeler & Meier 2002). 
The process of field generation and expulsion will depend on the details of the differential rotation 
(Akiyama et al. 2002; Akiyama, Wheeler & Lichtenstadt 2002). 

For this “wrapping” mechanism, the field will grow linearly with time and be proportional to 
the seed field. After revolutions of the neutron star, the initial seed (poloidal) field will produce 
a toroidal field, 

~ 2Tm^Bp, ( 8 ) 

where Bp is the initial seed poloidal field. The field will be limited by buoyancy, which will operate 
to expel the field from the site where it is generated. The field at this buoyancy limit, B^, is given 
by, 

h 2 


( 9 ) 
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where fb is the fractional difference in density between the rising flux tube elements and the stellar 
material. 

There are two regions where there is likely to be large shear and hence field growth. One 
is within the neutron star. The angular velocity can be large there, but, as remarked above, the 
gradient is rather flat after core bounce. A complication is that the dynamics of the LW phase 
may already have amplified the field within the neutron star so the seed field for the subsequent 
wrapping phase may be larger. The other region where strong field amplification is expected is at 
the boundary of the protoneutron star where the angular velocity is less than in the neutron star, 
but the difference in the angular velocity with respect to the postshock, settling matter results in 
a rather large gradient. 

After bounce, the central density of the neutron star is about 2 x 10^^ g cm The boundary 
of the homologous core at R ~ 10 km is initially at about 10^^ g cm at bounce, but settles to 
somewhat higher densities, ~ 3 x 10^^ g cm after some relaxation (Burrows, Hayes & Fryxell 
1995; Janka & Muller, 1996; Mezzacappa et al. 1998). For sound speed c* ~ lO^'^ cm s“^, density 
p ~ 10^^ g cm and fb ~ 0.01 the critical value of the field at the onset of buoyancy is, 

Bb-Qx 10^6 G (10) 

where fb -2 = /fe/0.01 and pu = p /lO^^ g cm The number of revolutions to reach the buoyancy 
limit Bb is, 

^ ^ ( 11 ) 

For Bq ~ 10^^ G and p r\j 10 ^^ g cm the linear amplification time scale before the field is 
expelled by buoyancy is thus, 

tf ~ UfP ~ 10 s Pms II-IpIa, (12) 

where P^s is the rotational period in msec. 

The second region, around the boundary of the protoneutron star at R ~ 50 km, has a typical 
density ~ 10^^ g cm after the shock stalls at about 200 km. For sound speed ~ 3 x 10® cm s“^, 
density p ~ 10^^ g cm and fb ~ 0.01, the critical value of the field at the onset of buoyancy at 
the protoneutron star boundary is, 

Bb^5x lO'" G flLlp\{\ (13) 

the number of revolutions to reach the buoyancy limit is, 

Uf 80 B^}, fl!_lpl{\ (14) 

and the linear amplification time scale before the field is expelled by buoyancy is, 

tf ~ 0.08 s Pms Pop 2 f&n^- 


(15) 
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Linear wrapping, therefore, is slow at first, but is likely to accelerate with the contraction and 
spin up of the neutron star. Because there will be a gradient in density and sound speed from within 
the neutron star to the post-shock, settling material, there will be variations in the time at which 
the field reaches the buoyancy limit and the value of the field at that limit around the surface of the 
neutron star. The magnetic field will be subject to wrapping by differential rotation throughout the 
collapsing medium since the field is ultimately anchored in the outer, slowly rotating progenitor. 

In general, the epoch of maximum expulsion of the toroidal magnetic field is seconds for the 
strong field within the neutron star to a fraction of a second somewhat beyond the neutron star, a 
range that spans the cooling, contraction time of ~ 1 s. An important point is that this process of 
wrapping and expelling field does not necessarily halt after one characteristic time. The toroidal 
field will build up to the buoyancy limit and then begin to be expelled. Field of this strength, 
^ 10^® G near the homologous core boundary, ^ 10^^ G at the boundary of the protoneutron star, 
will continue to pump out as the differential rotational energy is dissipated with a time scale that 
is much greater than the characteristic time of the differential rotation ~ 11“^. 

This discussion illustrates that the buoyancy time and subsequent, longer term dissipation of 
the rotational energy into magnetic field can be short enough to be a significant effect on the radial 
contraction time of the neutron star. There will be a gradient of buoyancy field that breaks out of 
the neutron star as a function of time and location, the nature of which needs to be explored. The 
characteristic timescales for growth and buoyancy throughout the neutron star are long compared to 
dynamical times and hence all the associated MHD processes will be quasi-static on that timescale. 
The energy in the magnetic field therefore does not directly power the explosion; that is, the 
magnetic pressure does not build up to explosive values and unbind the star. Rather, the magnetic 
field serves as a conduit to transmit the energy of differential rotation into directed jet energy 
and momentum over many dynamical times, and it is the jet dynamical power that creates the 
explosion. 

The rotational amplification of the field and its expulsion will go in the direction of eliminating 
the shear, but the infall of more slowly rotating matter will always maintain a differential rotation 
with negative gradient in fl at the neutron star boundary. The creation and expulsion of field from 
this layer should continue until the explosion of the star cuts off the infall. Initially, the magnetic 
field should have negligible feedback on the dynamics; but, as the magnetic field builds up, it will 
affect the inflow. 


4-2.2. Spring and Fling in Supernovae 

We can now assess the affect of the toroidal field on the physics of jet formation in the supernova 
ambience. The Blandford-Payne power (eqn. [2]) can be expressed as: 
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for parameters representative of the protoneutron star interior. Note that the effect of the smaller 
field near the protoneutron star boundary, ~ 10^^ G, would be substantially offset by the larger 
radius, ~ 50 km, so that the power is roughly the same in both key locations. If the neutron star 
can sustain a breakout buoyant field of 10^® G from deep within, or a field of 10^® G near the 
protoneutron star boundary, for ^ 0.1 s as the differential rotational energy is tapped, a strong 
supernova could result. 

For the conditions of interest here, the critical parameter of the magnetic switch mechanism 
(eqn. [5]) can be written as; 

1 / ~ 0.05 Si 5 (pnPmsf (17) 


where for convenience here and below, the mass of the neutron star has been normalized to 1.5 Mq. 
This equation says that there is a locus of critical values at which the switch would be activated 
and an especially high speed jet created, defined by i^=l, namely; 

B,r^t ^ 2 X 10^^^ G {puPmsf (18) 

The value u = 1 corresponds to a (critical) MHD power of (eqn. [4]), 

Lcrit ^ 1-2 X 10^"^ ergs"Vii/An^6'^^^^5- (19) 


The conditions to make this switch in jet character could be found in collapse conditions for espe¬ 
cially rapid rotation and hence large field strength or if a strong field becomes buoyant and “floats” 
into lower density conditions, or as the inflow density declines. Note the potentially important role 
of the collimation of the jet as reflected in the parameter /ah- This parameter would be of order 
unity for a slow Blandford-Payne jet with a broad base and size of order the diameter of the neutron 
star in which case a rather large luminosity is required to trigger the fast jet. If, however, the jet is 
tightly collimated as numerical simulations show for the case with “switch on” or as indicated by 
7 -ray burst jet collimation, then /ao could be quite small ^ 10“^ and a significantly more modest 
luminosity would be adequate to trigger the switch and the fast jet. Hoop stresses in the field 
surrounding the neutron star could lead to such tight collimation even for the slow jet, but this 
requires more study. 


The speed of the jet when the switch is “on” will be somewhere between the Alfven speed 

2 \ 1/2 
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and the “equipartition” speed (where Lmhd — Lkinetic) 

VE ^ ~ 1.7 X 10® cm s~^p^^^^p-^^ , 


( 20 ) 


( 21 ) 


which is ~ VA if, as we expect, the field amplification saturates and va ~ Hi?. These numbers 
suggest the possibility of mildly relativistic flows if B exceeds 10^® G, with the outflow speeds 
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depending critically on the maximum Alfven speed attained in the region outside the rapidly 
rotating core. These aspects clearly need further study. 

These order of magnitude estimates serve to illustrate that this spring mechanism might be 
quite common in core collapse, generating bipolar outflows of sufficient energy and velocity that 
they could be the sole means of supernova explosion. At minimum, the presence of such a jet would 
distort an already successful explosion, thereby creating the observed asymmetries. 


4-2.3. Dynamics and Nucleosynthesis 

There is a danger that the process we have outlined here is too effective in one regard, the 
production of neutron-rich elements. We have noted that this is a severe constraint on the LeBlanc- 
Wilson jet, such that the LW phase cannot be responsible for routine supernovae without severely 
over-producing neutron-rich species (§4.1). If conditions are just right, a jet-like explosion might 
account for the r-process (Symbalisty, Schramm & Wilson 1985; Cameron 2001). For an explosion 
based on the protopulsar jets, we have to seek a balance between an effective deposition of energy, 
namely a jet that is not too fast or narrow, and one that is not too broad that will lift off too much 
neutron-rich matter from deep in the collapse environment. Resolution of this dilemma will require 
more quantitative analysis, but here we sketch some of the basic issues. 

There is a major advantage to a jet that forms near or beyond the protoneutron star boundary 
where the neutronization will not be as extreme as it is deep within the protoneutron star where 
the LW jet forms and where large toroidal fields will form in the subsequent wrapping, contraction 
phase. Nevertheless, the infalling matter will still be neutron rich (the electron fraction Yg plum¬ 
mets from ^ 0.5 at the standing shock to ^ 0.2 for R ^ 50 km; e.g. Mezzacappa et al. 1998), 
and excessive ejection would again violate limits on neutron-rich species. The broad fan of the 
Blandford-Payne MHD wind that is conducive to a robust explosion thus may also be counter¬ 
productive from a nucleosynthesis point of view. There are several factors that may mitigate this 
production of neutron-rich ejecta. 

Hoop stresses and the large ambient pressure may restrict the opening angle of the initial jet. 
This would decrease the mass of neutron-rich matter directly ejected in the jet. We note that the jet 
should broaden as it passes though the standing shock at about 200 km due to the sudden decrease 
of the external pressure. The external pressure drops by about an order of magnitude at that 
point whereas the pressure of the plasma confined within the jet remains the same. The opening 
angle of the jet might then expand by about the ratio of the sound speeds, or a factor of about 3. 
This may help to broaden the fan of the wind, but at larger radii where matter is somewhat less 
neutronized. On the other hand, a broad jet fan may encompass all the very neutron-rich matter 
originally expelled in an LW bubble leading to expulsion of all that neutron-rich matter even if the 
LW phase did not result in an explosion. 

Another factor is that the jet is likely to go up the rotation axis where the density will be 
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relatively low due to the effects of rotation. After bounce, rotating models give densities at 100 
km about a factor of 5 - 10 lower on the rotation axis compared to the equator (e.g. Symbalisty 
1984; Monchmeyer, Schafer, Muller & Kates 1991). That means for a given opening angle of the 
slow toroidal jet, there will be relatively less matter to be intercepted by about the same factor. 
The attendant lower pressure will also tend to promote a broader fan, so this must be checked 
quantitatively. The rotation will also hinder convection in the equatorial direction and this may 
promote neutrino emission up the rotation axis (Fryer & Heger 2000). This tendency may also help 
to decrease the density along the rotation axis, the direction of propagation of the jet. 

We also note that the complementary effect is a higher density along the equator. For extreme 
values of rotation, this matter will form an accretion torus with an accretion time long compared to 
the dynamical timescale. The large rotation and distorted geometry of this protopulsar phase will 
evolve through shedding of angular momentum to the jet and environment to the more familiar, 
more nearly spherical, pulsar phase. 

Further out in the star, the rotation will cause less of an effect on the isodensity contours, but 
there will still be some distortion. In spherical core collapse supernovae models, iron is produced 
and ejected by the action of the shock on the silicon layers and the associated production of ®®Ni 
that subsequently decays to form ®®Fe. After several seconds, this silicon will be part of the infalling 
matter. The jet will create bow shocks that converge on the equator and drive overlying matter 
out preferentially in the equatorial plane (Khokhlov et al. 1999). The bow shocks will thus tend to 
drive the silicon layer and the newly formed nickel outward along the equator. This effect would 
be enhanced by any rotational distortion. Some of the pre-existing silicon may be directly ejected 
in this way and some silicon will be produced and ejected by shocking the oxygen layer. Silicon 
produced in the explosion may thus also be blown off in an equatorially preferential manner, but 
at somewhat lower velocities. These dynamics may partially account for the observation of high 
velocity iron lying outside lower velocity silicon within the body of Cas A (Hughes et al. 2000; 
Hwang et al. 2000) and for high velocity knots of silicon found beyond the reverse shock (Fesen 
& Gunderson 1996). The jets themselves should be formed of neutron-rich iron-peak material, 
but how much of this matter escapes the star will depend on the envelope structure and shock 
dynamics. 


4.3. The Magnetic Switch, Failed Supernovae, and the Gamma-Ray Burst Jet 

An important aspect of the dynamics of the jet propagation through a star is that the faster 
the jet, the less likely it is to trigger an explosion. According to the calculations of Khokhlov & 
Hoflich (2001) and Hoflich, Khokhlov & Wang (2001) (see also MacFadyen, Woosley & Heger 2001), 
faster jets have narrower opening angles, allowing less time for transverse shocks to affect the star 
before the jet emerges from the boundary of the core. The jet energy will tend to stay focused and 
directed along the jet trajectory rather than powering the bow shock and its associated transverse 
sheath of shocks that are critical for spreading the jet energy throughout the mantle and exploding 
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the star. (The difference is akin to the result of aiming a baseball and a needle with the same 
energy at a loaf of bread.) This leads to a suggestion for a different ultimate evolution for those 
conditions that lead to especially fast jets, namely a failed supernova, collapse to form a black hole, 
and the attendant possibility of creating a 7 -ray burst. These ideas are summarized schematically 
in Figure 2. 

We do not fully understand the conditions of either magnetic field strength or rotation of 
the progenitor core. It is plausible that these properties have a finite distribution such that a 
small fraction of progenitors have particularly large values of magnetic field and rotation. For a 
sub-set of all the core collapse events, e.g. in a given mass range, some could have properties 
that satisfy the criterion that the MHD luminosity exceeds the critical luminosity of the magnetic 
switch, Lmhd > Lcrit (§4.2.2). In this case, we expect the generation of an especially fast, narrow, 
jet with Vjet ~ VA, which may be well above Vesc- Such a jet may be less effective in creating a 
jet-induced supernova explosion. If this is the case, the mantle will continue to collapse and a black 
hole would form. The extreme version of this suggestion is that the narrow jet fails to form any 
supernova. It is plausible, however, that some matter is expelled to trigger a supernovae, but an 
insufficient fraction of the matter is ejected to prevent the ultimate fallback of enough mass to form 
a black hole. 

The delayed formation of a black hole would give a new ambience for the generation of a third, 
very possibly highly relativistic, jet, as widely discussed in the literature (MacFadyen & Woosley 
1999; Aloy, et ah 2000; MacFadyen, Woosley & Heger 2001; Zhang & Woosley 2001). Whether this 
relativistic jet can induce a supernova explosion when the earlier, slower jet failed to do so is an 
open question, but as we remarked above, the first jet could have generated a supernova through 
partial ejection of the mantle and still have allowed the formation of a black hole. In this case, the 
highly-relativistic jet would form within an already exploding configuration, an assumption made 
in some of the computations (e.g. Aloy et ah). 

What does seem clear is that if a new, relativistic jet can propagate out of the infalling matter, 
it will collide with the iron-rich matter expelled in the earlier jet. Depending on the timescale for 
the fallback to trigger the collapse to form the black hole (MacFadyen, Woosley & Heger 2000), 
the initial jet could have propagated well away from the star. For a typical velocity of the fast, 
but mildly-relativistic protopulsar jet of u ~ lO^'^ cm s“^, the matter will have reached a distance 
of order D ~ 10^^ ten ,2 cm, where tBH ,2 is the delay time between the initial collapse and the 
formation of the black hole and relativistic jet in units of 100 s. The material of the protopulsar 
jet may thus act as the “target” for the collision of the second jet and this collision could create 
the fireball and the internal shocks that represent the onset of the 7 -ray burst. The interaction 
between the 7 -ray burst jet and the target would occur at a distance D c/(c — u), which is of order 
D, but could be significantly greater as u —> c and the 7 -ray burst jet takes longer to catch up to 
the fast protopulsar jet. 

In the observer frame, the duration of the 7 -ray burst would be associated with the time for the 
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relativistic 7 -ray burst jet to propagate along the length of the mildly-relativistic jet, and this would 
be determined by the longer of two time scales; the duration of the non-relativistic, protopulsar jet 
(since this determines the length of the iron-rich lobe that is targeted) and the duration of the 7 -ray 
burst jet itself. The former time scale will be approximately the spindown time of the jet-emitting 
protopulsar: 

Tspindown — Ej~ot,NS/HD — 8 S Pms ^NS BnS,6' 

In this scenario where the tripping of the magnetic switch triggers an especially fast, inefficient jet, 
the field must exceed the critical value (eqn. [18]), so the spindown time will be 
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This spindown time can be as long as a typical 7 -ray burst, ~ 30 s, for rotation at less than breakup, 
tight collimation, and appropriate choice of density, so this could be the controlling factor in the 
duration of a 7 -ray burst if this scenario of colliding jets proves to have merit. 


Using a strong magnetic field as the prime carrier of momentum and energy in the second, 
relativistic jet, with the jet material just along for the ride, might explain how high Lorentz factors 
can be achieved and sustained in the black-hole-collapse jet. For example, for a “tenuous” corona 
around the accreting black hole of 10® g cm“^ in a magnetosphere of 10^® G, equations (20) and 
( 21 ) give relativistic Alfven/jet speeds of order 

/ n2 \ 1/2 

" ( 47 ^) == ““ 

The corresponding power from the Blandford-Znajek mechanism from a rapidly-rotating hole 
(Blandford & Znajek 1977; Meier 2001) is: 


Lgrb ^ 10 ®^ erg s ^ 


(25) 


where e is an efficiency factor and Mbh, 0.5 is the black hole mass in units of 3 Mq. If this luminosity 
is pumped into the relativistic jet for a time Trei seconds, the total energy emitted in the jet is; 

Eqrb ^ Lgrb Trei ^ 10 ®^ erg e Pel- (26) 

Recent analyses suggest that the typical energy emitted in 7 -ray bursts is around 3 x 10®® ergs 
(Kumar & Panaitescu 2001; Frail et al. 2001). If this low energy is typical, then very little of the 
expected black hole binding energy or rotation energy can be tapped to make 7 -ray bursts. For 
the picture we sketch here and, indeed, for any picture based on the production of jets from black 
holes we must have eB^^Trei ~ 0.03. In our particular magnetic switch scenario, we must satisfy 
equation (18) for the critical switch field and equation (24) to produce a suitably relativistic jet. 
This suggests that either the efficiency for the Blandford-Znajek process is relatively low (perhaps 
because the field is swallowed by the black hole rather than producing electromagnetic flux) or the 
duration of the process is short. We note that if the latter is the case, then the duration of the first 
jet (eqn. [23]) will control the duration of the 7 -ray burst. 
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Using the slower jet as the target and tying the delayed interaction between the two jets to 
the mantle collapse time would explain how 7 -rays are produced so efficiently and routinely far 
from the central black hole environment. In addition, the iron-rich nature of the first jet may give 
appropriate conditions in which to address the iron lines observed in some 7 -ray burst events (e. g. 
Piro, et al. 2000). The strength of the iron emission is difficult to assess. It depends on numerous 
very uncertain factors such as the mass of the iron ejected, the density in the precursor jet, the 
cross sections of the first and second jets, and the ambient ionizing flux. Estimates in the literature 
of the amount of iron required to produce a given line flux by recombination vary by orders of 
magnitude (e.g. Piro et al. 2000; Rees & Meszaros 2000; Meszaros & Rees 2001). 

This scenario of jet collisions does not help to resolve the issue of the low afterglow densities 
associated with some 7 -ray bursts (Kumar & Panaitescu 2001; Scalo & Wheeler 2002b) that are 
difficult to achieve in any model based on a massive star that expels a substantial wind. It also 
offers no explanation as to why the blast energy associated with the 7 -ray burst is confined to a 
narrow and rather small value, around 3 x 10^® ergs. The question of why such a small fraction of 
the available energy would be tapped in the jet is an important problem for the future. 


5. Summary and Discussion 

We have presented the case that strong toroidal magnetic fields will be generated in stellar 
collapse and that these magnetic fields can generate magneto-centrifugal jets in analogy to those 
found in simulations of black hole accretion. The case for magneto-centrifugal jets suggests they 
may be frequent and robust and hence provide a good basic understanding of why all core collapse 
supernovae are found to be substantially asymmetric and predominantly bi-polar. 

There are concerns that the jets we describe arise deep in the gravitational potential where the 
infalling matter is very neutron rich. We argue that this constraint rules out the prompt LeBlanc- 
Wilson jets as the common origin of supernovae. This suggests that conditions of initial rotation 
and magnetic field extreme enough to produce strong jets in this phase do not arise frequently 
in nature. We point out that successful jets from the later, protopulsar phase must walk the line 
between being broad enough to provide ample energy to the overlying matter without being so 
broad as to eject excessive neutron-rich matter. Partial evacuation of the matter on the axis by 
rotation may alleviate this problem, but it remains a concern. The jets will produce bow shocks 
that tend to expel matter, including iron and silicon, in equatorial tori. This may help to account 
for observations of the element distribution in Cas A. 

There also is some concern that gravitational radiation may remove so much angular momen¬ 
tum, and so quickly, that there would not be enough rotational energy to power a substantial MHD 
jet. We estimate that the effects of gravitational radiation (GR) will be significant, but not se¬ 
vere. The most important GR mode is probably the r-mode. Detailed computations by Lindblom, 
Tohline, and Vallisneri (2001) on the evolution and effects of this mode estimate that only about 



- 20 - 


half of the rotational kinetic energy will be lost before the gravitational radiation process switches 
off. The time scale for this to occur may be fairly short (a few tens of initial rotational periods, i.e. 
as short as, say, 30 ms), but the protoneutron star still will be left with a significant fraction of its 
initial rotational energy, which can be tapped for the production of MHD jets. 

The topics of core-collapse supernovae and 7 -ray bursts overlapped, at least in principle, with 
the discovery of SN 1998bw. Although the association is still controversial, this odd supernova is 
likely to have been connected to GRB 980425 (Galama et al. 1998). Supernova-like excesses of light 
have been detected in the afterglow of two 7 -ray bursts, GRB 970228 (Reichart 1999, Galama et 
al. 2000) and GRB 980326 (Bloom et al. 1999) about two weeks after the 7 -ray bursts. The excess 
light in the afterglow of both GRB 970228 and GRB 980326 has been modeled by the addition of 
light from an event like SN 1998bw. While there may be other explanations for these excesses, the 
connection with supernovae must be pursued. We suggest here a new way to make this connection. 
The conditions that apply in stellar collapse suggest that a magnetic “switch” mechanism found to 
apply in the black hole simulations may also apply in the collapse case with subsequent affect on 
the speed of the jet. If the switch turns on at low density and large magnetic field, an especially 
fast jet, with propagation speed of order the Alfven speed, could be produced. This jet would tend 
to propagate rapidly through the star and deposit relatively little momentum. The result might 
be a supernova explosion, but with enough infall to produce a delayed black hole. A reprise of 
the physics outlined here, and that explicitly invoked to produce relativistic MHD jets from stellar 
mass and AGN black holes, could then produce a relativistic jet that catches up to the first jet 
after it has emerged from the star. The interaction of these two jets could plausibly be the origin 
of the internal shocks thought to produce 7 -ray bursts and could explain the presence of iron lines 
in the afterglow. 

In this paper we have discussed toroidal fields of order 10^^ to 10^® G arising routinely in 
the formation of a neutron star in a supernova. There is clearly an issue of the final effective 
dipole strength of the neutron stars left behind in the explosion. We must have a final dipole field 
distribution of the neutron stars consistent with those of observed pulsars, ~ 10^2 to 10^3 G. It 
seems plausible that the large fields generated in the collapse are dissipated in the helical jets. The 
very strong fields that may be generated in the depths of the neutron star will take longer to float 
out, but the time still could be short compared to the lifetime of the supernovae, never mind the 
pulsar. Stronger fields in the interior could also leave smaller surface dipole fields. Much of the field 
may be amplified and dissipated at the shearing boundary between the infalling matter and the 
neutron star. This layer and its field will almost surely substantially dissipate after the explosion. 
These issues should be given more careful consideration. 

All of the issues raised here should be examined in the context of magnetars. However the large 
effective dipole fields arise in magnetars, there is a strong presumption that they are created in the 
formation of the neutron star. The formation of a magnetar may require special circumstances, 
or they may arise as the tail of the normal pulsar formation process. One possibility is that they 
represent that fraction of core-collapse events that, perhaps by accident of progenitor conditions, do 
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attain an especially rapid differential rotation and amplify the field not just by field line wrapping, 
but by an a —17 dynamo (Duncan & Thompson 1992). This could lead to exponential field growth, 
and total field strength of order 10^^ G with dipole component of order 10^^ G. This is approximately 
the regime of fast toroidal jets that we invoke to produce less ejecta and more infall to produce a 
black hole and perhaps a y-ray burst. To reconcile these pictures, it may be that magnetars arise 
with somewhat higher than normal fields and rotation, but still in the slow toroidal jet phase with 
conditions where the magnetic switch is still “off.” Even more extreme conditions could lead to the 
flipping of the switch, producing an especially fast jet that would be less effective in driving the 
explosion and thus leaving a black hole. Clearly there is a great need to study the many aspects of 
this problem in numerical and quantitative detail. 

Some models of SN 1998bw invoked especially large kinetic energies, in excess of 10^^ ergs in 
spherically-symmetric models, to account for the bright light curve and high velocities (Iwamoto, et 
al. 1998; Woosley, Eastman, & Schmidt 1998), while others took note of the measured polarization 
to suggest that strongly asymmetric models could account for the observations with more “normal” 
energies (Hoflich, Wheeler, & Wang 1999). Large kinetic energy has also been attributed to several 
other supernovae (again based on spherically symmetric models), especially SN 1997cy (Germany et 
al. 2000) and SN 1997ef (Branch 2001; Iwamoto et al. 2000). Events like SN 1997cy, SN 1997ef and 
SN 1998bw will help to sort out the physics of explosive events, whether such events are more closely 
related to “ordinary” supernovae or “hypernovae”, whether either of these classes leaves behind 
neutron stars as “ordinary” pulsars or highly magnetized “magnetars” or whether the remnant is 
a black hole and whether any of these events are associated with classic cosmic y-ray bursts as 
suggested by the supernova-like brightening of the afterglow of GRB 970228 and GRB 980326. In 
terms of our magnetic switch mechanism, it may be that “hypernovae” are still in the switch “off” 
phase so they have strong jets, but not so fast and narrow that the explosion is mitigated. In this 
sense, there may be a connection between the “hypernovae” and magnetars that, as argued above, 
may arise in the same regime of especially energetic protopulsar jets, but still with the switch “off.” 

If it were not for the constraints of the supernova polarimetry, we might continue to study 
essentially spherically symmetric collapse models, as indeed, many people will and must. It may be 
that we are pessimistic and that neutrino asymmetries can account for the observed polarimetry. 
The fact that jets alone can account for supernova explosions and the observed asymmetries has 
been established in principle. We think nature is saying that jets must occur in routine supernovae. 
We also believe that the process of neutron star formation must inevitably involve the physics 
outlined here whether the resulting jets are weak or strong. Glearly, we need a more rigorous 
description of the process of field generation, buoyancy, emergence, and subsequent evolution. 
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Fig. 1.— Schematic illustration of the protopulsar or slow toroidal jet phase. The left panel (a 
blow-up of the very inner portion of the right panel) shows the weak LeBlanc/Wilson jet emerging 
from deep within the newly born neutron star. Also shown is the torus of magnetic field that 
will be created by wrapping of the initial field that is anchored in the outer core. The wrapping 
leads to held amplihcation, buoyancy, and emergence of an accelerating helix of held (solid lines) 
that expands, carrying plasma in jet-like how driven by held torsion and conhned by held tension. 
The right panel shows the effect of the jet on the surrounding star. Bipolar jets propagate up the 
rotation axis. Bow shocks lead to the lateral propagation of shocks that collide on the equator and 
drive preferential motion outward there. The net effect is to drive a strongly asymmetric explosion 
into the outer portions of the star. 

Fig. 2.— Schematic illustration of the effects of a fast toroidal jet. Left panel (a blow-up of the very 
inner portion of the central panel) shows a tightly collimated, narrow, fast jet that would occur 
for conditions when the “magnetic switch” was “on” (see text). The fast jet can break out of the 
star, so most of the energy is dissipated along the jet, rather than in lateral shocks (see Figure 1). 
The result can be the failure to explode the star or a weakened explosion that results in continued 
collapse to form a black hole. The “failed supernova jet” will proceed to propagate outward, and 
the extreme conditions of rotation and gravity in the environment of the black hole can create an 
especially fast, highly-relativistic MHD jet (right panel). The relativistic jet can collide with the 
previously expelled matter of the “failed supernova jet” which acts as a target to create internal 
shocks and to produce a y-ray burst. 
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